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The  zero-one  mixed  integer  programming  problam  may  be  written 

Maximize  cx  +  dy 

aubject  to  Ax  +  Dy  <  b 

(1)  x  <  ef 

m 

x,  y  >  0  and  x  Integer 

where  A  la  mxn,  D  la  mxr,  x  la  nxl,  y  la  mxl,  and  c,  d,  b,  and  e  are  dimensioned 
compatibly.  Adding  vectors  u  and  v  of  alack  variables,  the  constraint*  of 

(1)  become 

Ax  +  Dy  +  u  -  bt+ 

(2)  x  +  v  -  e 

x,  y,  u,  v  >  0  and  x  integer. 

A  strategy  that  la  often  used  for  solving  (1)  is  to  adjoin  additional 
constraints  and  variables  to  partition  the  feasible  region  into  subsets  (e.g. , 
restricting  cx  +  dy  t<r  specified  Intervals),  and  seek  a  feasible  solution  to  the 
additionally  constrained  problem.  We  aosume  that  such  constraints  and  variables 
are  already  Included  in  (1)  and  (2),  and  address  ourselves  to  obtaining  feasible 
solutions  to  (2). 

In  a  variety  of  practical  situations  the  constralncs  of  (2)  contain  imbedded 
network  problems  (such  as  transportation  and  assignment  problems)  having  the 
property  that  every  extreme  point  solution  is  integer.  Moreover,  one  can  readily 
derive  constraints  from  (2)  that  impose  bounds  on  nested  partial  sums  of  variables 
(see  [3]),  and  these  constraints  also  have  the  property  that  every  extreme  point 
is  integer.  Finally,  a  set  of  such  constraints  can  be  adjoined  tc  further 
partition  the  feasible  region  of  (2). 

t  e  dei  es  a  vector  of  ones. 

tt  The  results  to  follow  are  unchanged  if  this  constraint  is  more  generally 
Ax  +  Dy  +  uU  -  b,  provided  the  augmented  matrix  (D,  U)  contains  an  mxm  non¬ 
singular  submatrix. 


2. 


Thus,  representing  these  special  constraints  by  Px  ■  £  and  Qx  +  w  -  g , 
consider  the  aupnented  system 

Ax  +  Dy  +  u  ■  b 


Qx  +  w  -  g 

x,  y,  u,  v,  w  ■  0  and  x  Integer, 

where  P,  Q,  £  and  g  are  Integer  matrices,  P  Is  pxn,  Q  Is  qxn,  and 


(0 


has  the 


unlmodular  property;  l.e.  ,  every  square  submatrix  of^^jhas  determinant  0,  1, 
or  -1.  (P  or  Q  may  also  be  null.  )  We  also  stipulate  p  <  n  and  every  pxp 
submatrix  of  P  Is  nonsingular. 

Theorem.  If  there  Is  a  feasible  solution  to  (3)  with  x  *  x'  and  x'  Integer, 
then  there  Is  a  basic  feasible  solution  to  (3)  with  x  »  x'  and  m  of  the 
components  of  (y,  u)  basic.  Moreover,  every  basic  feasible  solution  to 
(3)  with  m  of  the  components  of  (y,  u)  basic  assigns  Integer  values  to 
the  components  of  x,  v,  and  w. 

The  chief  significance  of  this  theorem  Is  that  it  permits  one  to  elect  a 
strategy  for  solving  (2)  that  focuses  on  finding  a  solution  to  (3)  with  m  of 
the  components  of  (y,  u)  basic. 

The  first  application  of  such  a  strategy  occurs  in  the  pure  zero-one  linear 
programming  method  of  Cabot  and  Hurter  [1],  whose  ideas  motivate  this  note. 
Specifically,  the  Cabot  and  Hurter  method  results  (for  D  and  d  null)  by 
adjoining  the  constraint  ex  -  N  to  (2)  and  replacing  b  by  b  +  ee,  where  0  <  e  <  1 
and  A  and  b  are  assumed  integer.  Beginning  with  N  at  an  upper  bound  for  ex, 

Cabot  and  Hurter  prescribe  finding  a  basic  feasible  solution  to  this  particular 
version  of  (3)  with  all  components  of  u  basic.  (No  procedure  is  given  for 


3. 


accomplishing  this,  however.)  If  en  acceptable  solution  is  found,  the  method 
stops  (or  is  applied  to  a  new  system  (2)).  Otherwise,  N  is  decremented  and  the 
process  repeats. 

The  theorem  implies  that  ex  ■  N  can  be  dispensed  with  in  the  Cabot  and 
Rurter  approach,  making  it  unnecessary  to  reapply  the  process  for  different 
values  of  N.  Moreover,  a  variety  of  other  side  conditions  and  supplementary 
constraints  can  be  accommodated  by  the  theorem,  some  of  the  more  important  of 
which  have  already  been  indicated. 

A  method  for  exploiting  the  theorem  for  the  pure  Integer  programming  problem 
is  given  below.  We  first  establish  the  validity  of  the  theorem  with  the 
following  three  lemmas. 

Lemma  1.  Every  square  submatrix  of 

(i  1  o\ 

H  ■  P  0  0  I  has  determinant  0,  1 ,  or  -  1 

°  v 

(where  the  I  matrices  of  the  top  row  are  nxn). 

Proof:  First,  it  is  assumed  that  has  the  unlmodular  property.  Using 
Induction  on  q,  and  expanding  the  determinants  of  the  appropriate  submatrices 
of  H  by  minors,  it  is  easy  to  see  that  ^  has  the  unlmodular  property.  The 
rest  of  the  proof  follows  the  same  argument,  using  induction  on  n. 

Lemma  2.  Consider  the  system 
Mt  +  Rx  -  ® 

Hz  -  p  and  t,  z  -  0 

where  M  is  mx-t,  R  is  mxs,  and  H  is  hxs,  the  vectors  a  and  f}  dimensioned 
compatibly.  If  f)  is  Integer  and  every  square  submatrix  of  H  has  determinant 
0,  1  or  -1,  then  z  is  integor  in  every  basic  solution  to  (4)  with  m  of  the 
components  of  t  basic. 
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Proof:  The  beeic  solution  must  have  the  form 

t  -  Mj*  (cr  -  and  z  -  H^p,  where  is  an  uocm  aubmatrix  of 

M,  Rj  la  an  mxh  aubmatrix  of  R,  and  la  an  hxh  aubmatrix  of  H.  The  latter 
assures  z  Is  Integer  (see  Hoffman  and  Kruskal  [5]). 

Lemmas  1  and  2  collectively  imply  the  latter  part  of  the  theorem.  The  first  part 
of  the  theorem  is  implied  by  the  following  stronger  statement. 


Lemma  3.  If  there  is  a  feasible  solution  to  (3)  with  x  ■  x'  and  x'  integer  then 

there  is  a  basic  feasible  solution  to  (3)  in  which 

(i)  every  component  of  w  is  basic, 

(ii)  Xj  is  basic  if  xj  ■  1  and  v^  is  basic  if  x!  ■  0 

(iii)  any  p  of  the  remaining  components  of  (x,v)  are  basic, 

(iv)  m  of  the  components  of  (y,  u)  are  basic. 

Proof:  Conditions  (i),  (ii),  (iii)  give  q  +  n  +  p  variables.  We  show  that  the 

submatrix  composed  of  the  associated  columns  of  H  of  Lemma  1  (call  it  H  )  is 

nonsingular.  First,  all  columns  of  Q  in  H,  may  be  reduced  to  0  by  subtracting 

from  them  appropriate  multiples  of  the  columns  of  1^  (associated  with  w). 

Next,  there  must  be  exactly  p  indices  j  such  that  x^  and  v^  are  chosen 

to  be  basic.  Subtracting  each  of  these  v  columns  from  its  associated  x  column, 

J  /p  R\  J 

and  rearranging  rows  and  columns,  transforms  into  (  1  J  ,  where  P^  is  a  '  J 
pxp  submatrix  of  P,  R  consists  of  0  columns  and  columns  of  P,  and  I  is  the 
(tH-q)x(nfq)  identity  matrix.  The  nonsingularity  of  asures  is  non- 

singular,  and  hence  a  basis  for  the  subsystem  of  (3)  with  Ax  +  Dy  +  u  •  b  removed. 

Since  contains  a  column  for  each  positive  component  of  x',  t',  and  w' 

(where  t'  ■  e  -  x* ,  w'  -  g  -  Qx'),  the  basic  solution  to  the  subsystem  must  yield 
x  -  x'.  Finally,  (iv)  is  established  in  conjunction  with  (i),  (ii),  and  (iii) 


from  the  assumed  existence  of  a  feasible,  and  hence  a  basic  feasible,  solution 
to  Dy  +  u  ■  b  -  Ax' . 


V 


. ..  .K.  rrtnw 


ur i'**' wi  /  f*  ■n  ffl* ^  W  V^.V^I? 


5. 


An  Integer  Programming  Method 

We  give  an  Integer  programming  method  for  the  pure  aero -one  problem  (with 
D  and  d  null)  that  pursues  the  objective  of  making  the  m  components  of  u  basic 
in  (3).  Assume  A  and  b  are  integer,  and  replace  A  by  2A  and  b  by  2b  4-  e. 
(This  replacement  clearly  does  not  change  the  set  of  integer  x  satisfying 
Ax  <  b,  and  implies  u  >  e  and  integer  for  all  non -negative  u,  x  satisfying 
u  +  Ax  ■  b  and  x  integer.  ) 

1.  Solve  the  linear  program 

Maximize  x  -  au 
o 

subject  to  (3)  (disregarding  the  integer  restriction  on  x), 
where  a  >  0  and  integer  (e. g. ,  let  a  «  e). 

Represent  the  current  tableau  for  the  simplex  method  in  the  form 


% 


Maximize  x  -  a  +  £  a  .  (-t.) 

0  °°  j-1  °J  J 


'1  •  ‘to +  £  -U  <-V  1  ■ 1 . - 

where  the  are  the  current  basic  variables  and  the  t^  the  current  nonbasic 

variables.  (The  a^  coefficients  of  the  current  tableau  are  not  to  be  con¬ 
fused  with  the  components  of  the  A  matrix. ) 

Upon  obtaining  an  optimal  tableau  (a^Q  >  0  and  a  ^  >  0  for  i,  J  >  1), 
go  to  Step  2. 


2.  If  a.  is  integer  for  all  1,  the  basic  solution  z.  -  a. 

io  i  io 

and  t^  ■  0  (all  i,  j)  gives  a  feasible  integer  solution  for  (2)  by  identifying 
the  variables  x^  from  among  those  currently  designated  z^  and  t^.  Otherwise, 
if  some  a^o  is  noninteger,  adjoin  a  cut  [4,  2]  and  re-optlmize  with  the  dual 
simplex  method. 


1 


>A«  *-». 
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3.  If  same  is  still  noninteger,  let  t^  denote  the  current 

nonbsslc  variable  that  was  the  slack  variable  for  the  cut  adjolnted  (most  recently) 
in  Step  2.  Replace  a  with  a  -  K  <  0,  where  K  is  an  Integer  (e.g.  ,  the 
least  Integer  <  Then  re-optimlze  with  the  primal  simplex  method  and 

return  to  2. 

The  purpose  of  Step  3  is  to  exploit  the  fact  that  the  cut  slack  qualifies 
to  be  one  of  the  "u^"  of  the  theorem.  Thus,  it  assigns  the  slack  a  weight  in 
the  objective  function  designed  to  drive  it  basic,  thereby  possibly  modifying, 
but  not  discarding,  the  weights  assigned  to  the  other  u^. 

Finite  convergence  is  guaranteed  if  one  uses  the  choice  rules  of  Gomory  [4] 
in  Step  2  and  bypasses  Step  3  after  a  fixed  number  of  iterations.  The  main  point, 
of  course,  is  that  the  method  gives  a  way  to  pursue  integer  feasibility  by 
exploiting  the  theorem. 
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